Snake states and Klein tunneling in a graphene Hall bar with a pn-junction 
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The Hall (Rh) and bend (Rb) resistance of a graphene Hall bar structure containing a pn-junction are 
calculated when in the ballistic regime. The simulations are done using the billiard model. Introducing a pn- 
j unction-dividing the Hall bar geometry in two regions-leads to two distinct regimes exhibiting very different 
physics: 1) both regions are of n-type and 2) one region is n-type and the other p-type. In regime (1) a 'Hall 
plateau'-an enhancement of the resistance-appears for Rh- On the other hand, in regime (2), we found a 
negative Rh-, which approaches zero for large B. The bend resistance is highly asymmetric in regime (2) and 
the resistance increases with increasing magnetic field B in one direction while it reduces to zero in the other 
direction. 



PACS numbers: 72.80.Vp, 73.23.Ad, 73.43.-f 



In graphene, a single flat layer of hexagonal carbon 
atoms, the carriers act as a 2D ultra-relativist ic fermion 
gas. Among interesting differences with standard two- 
dimensional electron gases (2DEGs) we should mention: 
the extra-ordinary Quantum Hall effect, Klein paradox 1 , 
and Zitterbewegung 2 , as some of the most striking char- 
acteristics of this new material. Although for most of 
these effects only indirect signatures were found in exper- 
iments up to now, graphene's behaviour as a relativistic 
analogue has far-reaching consequences in electronic de- 
vices when based on this material. In this context, the 
pn-junction is an interesting system because it exhibits 
Klein tunneling as well as the analogue of a negative re- 
fraction index. Moreover the phenomen of 'snake states' 
along the pn-interface has been predicted 3,4 . Actually 
at such a pn-junction one has a combination of 'skip- 
ping orbits' when electrons arriving at the pn-junction 
are reflected and 'snake states' when they exhibit neg- 
ative refractive index behaviour. Along with theoreti- 
cal proposals 4 , experiments on such systems were under- 
taken recently 5 ' 6 . This motivated us to investigate the re- 
sponse of a graphene Hall bar 7 ' 8 containing a pn-junction. 
In previous work on mesoscopic Hall bars Beenakker et 
al 9 proposed a classical model to describe the magnetic 
response. Such Hall probes have been applied as nonin- 
vasive probes for local magnetic fields 10 ' 11 to investigate 
e.g. mesoscopic superconducting disks and ferromagnetic 
particles. In this work we investigate a Hall-bar made of 
graphene with a pn-junction along one of its axes in the 
ballistic regime using the billiard model. 

The system we investigate, shown in Fig. 1, is a Hall 
bar with four identical leads, with in the middle a pn- 
junction dividing the Hall cross in an n-doped region and 
a p-doped region. The doping can be realized electrically 
by applying gate potentials V p and V n on the p- and 
n-region respectively. We modelled this by a shift in en- 
ergy Ep — > Ep — V n and define V = V p — V n: which 
is equivalent to have V n = 0. We use the following 





n 


p 

1 ^ _ 

\ 


\ 






/ 


1 
\ 






^1 

\ 
1 


a ^ -0 
/ 






1 

/ 

1 

) / ' 


\ 











V 



n 



FIG. 1. (Color online) A sketch of the Hall-bar geometry, 
where the pn-junction is represented by the red (dark grey) 
line, and the contacts are shown in blue (black) having a 
chemical potential \±i. Three configurations are possible for 
the voltage and current probes: a (Hall resistance measure- 
ment), and /3 and 7 (bend resistance measurements). The 
inset shows the potential profile. 



dimensionless units throughout the work, R* = R/Ro, 
with R = = B/Bq, with B = Jg^, 

r* = r c /W = (1 - 0/ B \ with f = V/E F , where W is 
the width of the channels, vp ~ 10 6 m/s is the Fermi- 
velocity, e is the charge of the electron. Our numerical 
simulations are done for typical values of W = 1/im and 
E F = 50meV, which results in R = 6.08 x 10 2 ^2 and 
Bq = 0.05T. In order for our approach to be valid, we 
must require that: (1) we are in the ballistic regime, 
meaning that the mean free path l e is larger than the 
width of the channels W, and (2) the system can be 
described classically, implying that, a) the confinement 
quantization E q = Tivp/W and b) the Landau level quan- 
tization Tiujd = \/2Tivf/Ib with the magnetic length 
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Ib = y/Ti/eB are small compared to the thermal fluc- 
tuations ksT. Since in graphene the temperature depen- 
dence of l e is found to be small within a broad range 
of temperature, and l e is of order of microns for clean 
graphene samples, see e.g. Ref. 8 , (1) is satisfied for not 
too low temperatures, and for the used parameters. 

We calculate the Hall and bend resistances according 
to the Landauer-Biittiker theory 13 where the resistances 
are defined as 
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with the voltage Vki = (pk — ^0/ e measured between 
terminals (fc, I) and a current I mn is flowing between 
(m, n). Tij are the transmission matrix elements giving 
the probability for an electron injected from terminal j 
to end up in terminal z, and D is a (sub) determinant 
of the transmission matrix T which is symmetric in the 
magnetic flux: D(B) = D{-B)). 

One has the following exact symmetries: (1) Rmn,ki = 
-Rmn,lk = -Rnrn,kh and (2) the reciprocity relation 
Rmn,ki(B) = Rki,mn(—B). There are three resistances 
of interest: the Hall resistance Rh = R a = ^13,24, the 
bend resistances Rb = Rp = ^14,32 and R 1 = i?i2,43, 
and their counterparts (denoted by a prime) obtained by 
switching the voltage and current probes, where the lat- 
ter are related to the former by the relation Ri(—B) = 
—R[(B). These resistances are defined by Ri — Vi/Ii 
(i = a,/3,7), as shown in Fig. 1. 

To obtain the transmission matrix elements we use the 
semi-classical billiard model 9 . are found by inject- 
ing a large number N of electrons from terminals j and 
tracking down which ratio ends up in terminals i. The 
probability to be injected from a terminal j is according 
to Ref. 9 uniformly distributed over the channel width, 
and has an angular distribution P(a) = ^cos(a). In 
graphene it also depends on the DOS which changes with 
the parameter £ = Ep/V . Because in the system under 
consideration the DOS is not the same in every chan- 
nel (and even not homogeneous inside channels 2 and 4) 
we incorporate this by a factor Sj(x) = |1 — £j(x)\ for 
injection from position x at terminal j. 

The trajectories of the electrons are assumed to be bal- 
listic and the electron moves along a circular orbit with 

cyclotron radius r c . Unlike in a standard 2DEG where 
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in graphene we have a linear dis- 
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persion relation E^ in = vfP and we obtain instead 
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a cyclotron radius which can be tuned by changing the 
magnetic field strength B or the potential V. To inves- 
tigate the behavior of an electron incident on the pn- 
junction, we calculate the transmission and reflection 
probabilities on the pn-j unction quantum mechanically 
(according to the Dirac model), as well as, how the par- 
ticle is deflected by the junction. The direction of the 



particle will change on transmission through the junc- 
tion. If the incident angle is fa then in case of reflection, 
the reflected angle is fa = fa. For transmission the angle 
fa of the transmitted electron is related to the incident 
one as 



p t sm(fa) = pi sin(fa), 



(3) 



where p t and Pi are the corresponding momenta in the 
respective regions, if the region has potential V then p = 
(E — V)/vp. In the case |sin(^)| > 1, total reflection 
occurs and the reflection probability R = 1. Otherwise 
R is given by 4 
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Thus when an electron hits the pn-junction, it has a prob- 
ability to be reflected or to be transmitted. 

Although the four-fold symmetry is lost in our system, 
we still have the following symmetries: 

1) Tij(B) — Tji(-B), the Onsager relation, 

2) T j2 (B) = T j4 (-B), with j = 1, 3, 

3) T 1S (B) = T 1S (-B), 

where the last two are specifically for our pn-junction 
set-up. Therefore, both the Hall and bend resistances 
obey 



1) R^-B) 

2) Ra(-B) 



-Rp{B), 

-Ra(B), 



Hence it suffices for the Hall resistance to only plot 
R a (B) for positive values of the magnetic field. More- 
over for the bend resistance we choose to plot Rp(B). 

To understand the influence of the pn-junction on the 
resistances we plot the electron density p(x, y) in Fig. 2 
for B* = 3, V/E F = 2. For this B the cyclotron radius 
r c = W/3 is equal in both regions, but with opposite di- 
rection of motion for the electrons. Notice that r c = W/3 
is large enough that the charge carriers will reach the pn- 
junction, but not the other edge of the channel without 
scattering (remember the width of the channels is W). 
There is symmetry for injection from lead 1 and 3 be- 
cause of equal cyclotron radius in both regions. For in- 
jection in the first lead we see that due to the negative 
refraction index there is a nonzero T21. Comparing in- 
jection from lead 2 and 4 we notice that the snake states 
are only travelling from 4 to 2 and not in the opposite 
direction. 

In Fig. 3 the (a) Hall resistance, (b) Hall coefficient, 
and (c) the bend resistance are shown for various val- 
ues of £ = V/Ep. Two cases are clearly distinct: (1) 
< £ < 1, such that the potential barrier is lower than 
the Fermi-energy and both regions are of n-type, and 
(2) £ > 1 where the right region is p-type, which re- 
verts the circling direction of the electrons and makes 
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FIG. 2. (Color online) The electron density p(x,y) inside the 
Hall cross for B* = 3, V/Ef — 2 for carrier injection from 
the four different leads as indicated by the arrows. 



the refraction index negative allowing for Klein tunnel- 
ing in graphene. From this figure we see that the Hall 
resistance shows for small < £ < 1, a plateau, while 
for £ > 1 it is negative and approaches zero with in- 
creasing £. The plateau in Rh is due to the pn-j unction 
which guides the electrons making it difficult for them to 
cross the junction 14 . If £ is increased above 1, the elec- 
trons in the p region will turn in the opposit direction, 
therefore they can bend towards a different channel and 
the pn-junction will collimate the electrons 15 diminish- 
ing the Hall resistance with increasing £ = V/Ep until 
it vanishes at V = 2Ep where r c is the same in both 
regions. From the Hall coefficient we can see that unlike 
the case without pn-junction, the Hall coefficient is neg- 
ative and decreases for high magnetic fields in the case 
£ > 1, while for £ < 1 it approaches asymptotically the 
value for £ = with increasing B. The bend resistance on 
the contrary becomes highly asymmetric, with increasing 
£ = V/Ep and is strongly enhanced by the presence of 
the pn-junction. In particular, notice that in regime (2) 
\Rb\ increases for large and negative B. To understand 
this we observe from Eq. (1) that Rb oc T31T24 — T34T21. 
Notice further that for large and negative £?, T31 & 
while T21 and T 34 are nonzero because the right region is 
of p-type. 

In summary, we investigated the Hall and bend re- 
sistance of a Hall cross made of graphene containing a 
pn-junction. The Hall resistance exhibits a 'last Hall 
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FIG. 3. (Color online) (a) Hall resistance, (b) Hall co- 
efficient, and (c) the bend resistance, versus the magnetic 
field strength are shown for various values of £ = V/Ef = 
0; 0.4; 0.8; 1.2; 1.6; 2, given by the labels next to the curves. 



plateau', as in a standard 2DEG when the Hall bar has 
rounded corners, in the case of £ < 1 where the junction 
serves as a guiding center for the electrons. For £ > 1, the 
Hall resistance is qualitatively very different: i) its sign 
is opposite, ii) no Hall plateau is present but rather a 
local minimum, hi) \Rh\ diminishes both with increasing 
£ and £?, and iv) for £ = 2, Rh = 0. The bend resistance 
is highly asymmetric in B for large V. 
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